IPPP/02/70 
DCPT/02/140 



On supersymmetric contributions to the CP 
asymmetry of the B —> <j)Ks 



m 

O ' S. Khalil 1 ' 2 and E. Kou 1 

O 

1 IPPP, Physics Department, Durham University, DH1 3LE, Durham, U. K. 

^ ' 2 Ain Shams University, Faculty of Science, Cairo, 11566, Egypt. 

>—> ' 

Abstract 

> 

£T) ' We analyse the CP asymmetry of the B — > (f)Ks process in general supersymmetric mod- 
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els. In the framework of the mass insertion approximation, we derive model independent 
limits for the mixing CP asymmetry. We show that chromomagnetic type of operator may 
play an important role in accounting for the deviation of the mixing CP asymmetry be- 
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tween B — > </>Kg and B — ► J/ipKs processes observed by Belle and BaBar experiments. 

A possible correlation between the direct and mixing CP asymmetry is also discussed. 

Finally, we apply our result in minimal supergravity model and supersymmetric models 

' with non-universal soft terms. 
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1 Introduction 

With the advent of experimental data from the B factories, the Standard Model (SM) 
will be subject to a very stringent test, with the potential for probing virtual effects from 
new physics. Measurements of the CP asymmetries in various processes are at the center 
of attentions since in the SM, all of them have to be consistently explained by a single 
parameter, the phase in the Cabbibo-Kobayashi-Maskawa mixing matrix 8km [!]• 

The BaBar [2] and Belle [3] measurements of time dependent asymmetry in B — ► 
J /ipKs have provided the first evidence for the CP violation in the B system. The world 
average of these results, Sj/^k s = sin 2/5(20!) = 0.734±0.054, is in a good agreement with 
the SM prediction. Therefore, one may have already concluded that the KM mechanism 
is the dominant source of the CP violation in B system. However, in the B — > J/ipKs 
process, new physics (NP) effects enter only at the one loop level while the SM contribu- 
tions has dominant tree level contributions. Thus, it is natural that NP does not show 
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up clearly in this process. In fact, in order that a significant supersymmetric contribu- 
tion appear to Sj/^Ksi we need a large flavour structure and/or large SUSY CP violating 
phases, which usually do not exist in most of the supersymmetric models such as SUSY 
models with minimal flavour violation or SUSY models with non-minimal flavour and 
hierarchical Yukawa couplings (see Ref. [4] more in detail). 

Unlike the B — > J/ipKs process, B — > <pKs is induced only at the one loop level both in 
SM and NP. Thus, it is tempting to expect that the SUSY contributions to this decay are 
more significant [5-8]. Based on the KM mechanism of CP violation, both CP asymmetries 
of B — > 4>Ks and B — > J/ipKs processes should measure sin 2/3 with negligible hadronic 
uncertainties (up to 0{\ 2 ) effects, with A being the Cabbibo mixing). However, the recent 
measurements by BaBar and Belle collaborations show a 2.7<x deviation from the observed 
value of Sj/^k s [3,9]. The average of these two measurements implies 

S+Kg = -0.39 ±0.41. (1) 

This difference between Sj/^k s an d S</>k s is considered as a hint for NP, in particular for 
supersymmetry. Several works in this respect are in the literature with detail discussion 
on the possible implications of this result [10-21]. 

As known, in supersymmetric models there are additional sources of flavour structures 
and CP violation with a strong correlation between them. Therefore, SUSY emerges as 
the natural candidate to solve the problem of the discrepancy between the CP asymme- 
tries Sj/fKs and S<pK s - However, the unsuccessful searches of the electric dipole moment 
(EDM) of electron, neutron, and mercury atom impose a stringent constraint on SUSY 
CP violating phases [22]. It was shown that the EDM can be naturally suppressed in 
SUSY models with small CP phases [22] or in SUSY models with flavour off-diagonal CP 
violation [22,23]. It is worth mentioning that the scenario of small CP phases ( < 10~ 2 ) 
in supersymmetric models is still allowed by the present experimental results [24]. In this 
class of models, the large flavour mixing is crucial to compensate for the smallness of the 
CP phases. 

The aim of this paper is to investigate, in a model independent way, the question 
of whether supersymmetry can significantly modify the CP asymmetry in the B — > <pKs 
process. We focus on the gluino contributions to the CP asymmetry S^,k s f° r the following 
two reasons. First, it is less constrained by the experimental results on the branching ratio 
of the inclusive transitions B — > X s j and B — > X s l + l~ than the chargino contributions 
[25]. Second, it includes the effect of the chromomagnetic operator which, as we will 
show, has a huge enhancement in SUSY models [26,27]. We perform this analysis at the 
NLO accuracy in QCD by using the results of Ali and Greub [28]. We also apply our 
result in minimal supergravity model where the soft SUSY breaking terms are universal 
and general SUSY models with non-universal soft terms and Yukawa couplings with large 
mixing. 

The paper is organised as follows. In section 2, we present the CP violation master 
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formulae in 5-system including the SUSY contribution. In section 3, we discuss the 
effective Hamiltonian for AB = 1 transition. Section 4 is devoted to the study of the 
supersymmetric contributions to the mixing and direct CP asymmetry S^x s and C$k s - 
We show that the chromomagnetic operator plays a crucial role in explaining the observed 
discrepancy between S^Ks and S^k s - In section 5, we analyse the SUSY contributions 
to S^Ks m explicit models. We show that only in SUSY models with non-universal soft 
breaking terms and large Yukawa mixing, one can get significant SUSY contributions to 
S(/>k s - Our conclusions are presented in section 6. 



2 The CP Violation in B -> <j)K s Process 

We start the sections by summarising our convention for the CP asymmetry in B system. 
The time dependent CP asymmetry for B — > (f)Ks can be described by [29]: 

r(B°(f) - <f>K s ) + T(B°(t) -> <f>K 8 ) 
= C^ Ks cos AM Bd t + S^Ks sin AM Bd t (3) 

where C^k s and S^k s represent the direct and the mixing CP asymmetry, respectively 
and they are given by 



The parameter ~p((f)Ks) is defined by 



(4) 



pms) = A^Kiy (5) 

where A(4>K S ) and A((f)K s ) are respectively the decay amplitudes of 5° and B° meson 
which can be written in terms of the matrix element of the AB = 1 transition as 

A( ( f>K s ) = ( ( j ) K s \nf B=1 \B^), A^K s ) = ^K s \nf B U\B"). (6) 

The mixing parameters p and q are defined by \B ± ) = p\B°) +q\lf), \B 2 ) = p\B°) — q\B°) 
where |-E>i(2)) are mass eigenstates of B meson. The ratio q/p can be written by using the 
off-diagonal element of the mass matrix and its non-identity q/p ^ 1 is the signature of 
the CP violation through mixing: 

q I Mf 2 - ir; 2 

P ^M 12 -|r 12 ' {) 
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The off-diagonal element of the mass matrix is given by the matrix element of the AB = 2 
transition as 

(B^nfsUB^^M^-'-Tn. (8) 

In SM, the major contribution to this matrix element is obtained from the box diagram 
with lU-gauge boson and top quark in the loop. As a result, we obtain: 

q = VttVtd (Q) 

v v tb v t * d { ) 

where we ignored terms 0(ri 2 /Mi 2 ). Since p((pKs) = ^ M j^| = v*v td = ^' ^ ae m i x i n S 
CP asymmetry in B — > 4>K S process is found to be 

S^ s = sin2/3. (10) 

Therefore, the mixing CP asymmetry in B — > <pKs is same as the one in B — > J/ifjKs 
process in SM. 

In supersymmetric theories, there are new contributions to the mixing parameters 
through other box diagrams with gluinos and charginos exchanges. These contributions 
to the AB = 2 transition are often parametrised by [30, 31] 



where M 12 = Mf^ 1 + M^ SY . In this case, the ratio of the mixing parameter q/p can be 
written as 

g - c -2i9 d V tb V td / 12 x 

p v tb v t * d y ' 

Thus, in the framework of SUSY, the mixing CP asymmetry in B — > J/^K S is modified 
as 

S jms = sm(2(3 + 26 d ). (13) 

In B — > 4>K S process, we have to additionally consider the SUSY contributions to the 
AB = 1 transition. The supersymmetric contributions to the AB = 1 transition comes 
from the penguin diagrams with gluinos and charginos in the loop (see Fig.l). We can 
parametrise this effect in the same manner [5,31]: 

A ^ Ks) S A e^, (14) 



A SM {(f)Ks) 



where A(<pK s ) = A SM (<pK s ) + A SUSY (<pK s ). Therefore, we obtain p(<pK s ) = e~ 2idA , 
hence Eq. (4) leads to 

<V S = 0, S 4>Ks = sm(2f3 + 29 d + 29 A ). (15) 
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However, this parametrisation is true only when we ignore the so-called strong phase. 
Since the Belle collaboration observed nonzero value for C^Ks [3] we should consider 
the strong phase in the analysis. In this respect, we reparametrise the SM and SUSY 
amplitudes as 

A SM {(t)Ks) = \A SM \e l5sM , A SUSY {<j)K s ) = \A SVSY \e idsuSY e lSsuSY , (16) 
A SM (<PK S ) = \A SM \e^M y A SUSY {<PK S ) = \A SVSY '| e -«W e iW ? (17) 

where Ssm(susy) is the strong phase (CP conserving) and Osusy is the CP violating phase. 
By using this parametrisation, Eq. (4) leads to 

sin 2(3 + 2 (^g^) cos S 12 sm{6 SUSY + 2(3) + (^g) ' sin{20 SU sY + 2(3) 

^ KS ~ \ c /MSUSYIX " ~ /| A SUSY|\2 

1 + 2 (^45^) C0S5 12 COSOSUSY + (^SMT 1 ) 

2 ( '^smI 1 ) sin S 12 sin 9 S usy 

<t> K S - ~ /|ASUSY|\ " " /| A SUSY|\2> V ) 

1 + 2 [ L T^w l ) cos 612 cos Osusy + [ ! ^swf ) 

where 5 12 = Ssm — 8susy- Assuming that the SUSY contribution to the amplitude 
is smaller than the SM one, we can simplify this formula by expanding it in terms of 
|A SUSY |/|A SM | [32]: 



I^susyj 
\A SM \ ' 

M SUSY 1 



^4>k s = sin 2(3 + 2 cos 2(3 sin S usy cos S 12 , /1SM| , (20) 



C^Ks = -2 sin Osusy sin 5 12 ^ SM ^ , (21) 

where (9((|A SUSY |/|A SM |) 2 ) is ignored. However, as can be seen by comparing these 
formulae to the measured value of sin 2(3 and the following Belle measurements 

S 4>Ks = -0.73 ±0.66, (22) 
C 4>Ks = -0.56 ±0.43, (23) 

a large value of |A |/|A SM | seems to be required. Therefore, in our analysis, we use 
the complete expressions for S^k s and C^k s as given in Eqs. (18) and (19), respectively. 

Finally, we comment on the branching ratio of the B° — ► (j)K° decay. Since we consider 
a reasonably large |A SUSY |/|A SM |, one may wonder that the prediction of the branching 
ratio could be significantly affected. The SUSY effect to the branching ratio can be 
expressed by 



Br(B° -> (f)K ) = Br SM (B° -> <f)K ) 



1 + 2 cos {Osusy - £12) -n; r + 



Asusy\ , [\Asusy\\ 



Asm\ V Asm 



(24) 
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where Btsm{B° — > <f)K°) is the standard model prediction of the branching ratio and 
obtained as (5 ~ 8) x 10 -6 . The branching ratio is measured by B factories as 

Br (B° -> <f)K Q ) = (10.0 ±\* x lO" 6 Belle [33] (25) 
Br(B° — > 4>K°) = (8.1+^ ±0.8) x 10~ 6 BaBar [34], (26) 

which in fact, imply that |A SUSY |/|A SM | can be of 0(1) (depending on the value of the 
phases Osusy and ^12) • As we will show in the following sections, this is the right magni- 
tude for enhancing S^k s - 



3 Effective Hamiltonian for AB = 1 transitions 

The Effective Hamiltonian for the AB = 1 transitions through penguin process in general 
can be expressed as 



n 



AB=1 

eff 



where 



G_f 
V2 



V tb V* 



CiOi + c g o g ]T cA + c g d 6 



li=3 



i=3 



(27) 



3 = S a YLb a SpYLSp, 

04 = Sa^LbpSp^LSa, 

5 = Sa^LbaSffLsp, 

06 = s a YLbps p ^Rs w 

_^£_ 

Ug - 8tt 2 



m b s a(J ^R-^bpG 



[IV 



(28) 
(29) 
(30) 
(31) 

(32) 



!5- 



The terms with tilde are obtained from Cj ifl and Oi t9 



where L = 1 — 75 and R = 1 
by exchanging L <-> R. The Wilson coefficient Cu g \ includes both SM and SUSY contri- 
butions. In our analysis, we neglect the effect of the operator 7 = -^rn h s OL a iil ' Rb^F^ 
and the electroweak penguin operators which give very small contributions. 

In this paper, we follow the work by Ali and Greub [28] (the generalised factorisation 
approach) and compute the B — > (pK process by including the NLL order precision for 
the Wilson coefficients 03^6 and at the LL order precision for C g . A problem of the gauge 
dependence and infrared singularities emerged in this approach was criticised in [35] 
and was discussed more in detail in [36]. Improved approaches are available these days, 
i.e. the perturbative QCD approach (pQCD) [37] and the QCD factorisation approach 
(BBNS) [38] and B — > <pK process has been calculated within those frameworks in 
[39,40] for pQCD and [41] for BBNS. However, in our point of view, the small theoretical 
errors given in [39-41] are too optimistic since both approaches are still in progress and 
for instance, the power corrections from either a s or l/m^ which are not included in the 
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theoretical errors of Refs. [39-41] could be sizable. Since the purpose of this paper is not 
to give any strict constraints on fundamental parameters of SM or SUSY, but to show if 
SUSY models have any chance to accommodate the observed deviation in CP asymmetry, 
we should stretch our theoretical uncertainties as much as possible. Therefore, we use Ali 
and Greub approach which include very large theoretical uncertainty in the prediction of 
the amplitude. In addition, we treat the CP conserving (strong) phase as an arbitrary 
value. Since a prediction of the strong phase is a very delicate matter, especially due to 
our ignorance of the final state interactions, we would like to be very conservative for its 
inclusion. 

The Wilson coefficient at a lower scale \i ~ 0(m&) can be extrapolated by 

Cii/j,) = U(n,n w )Ci(nw) i = l~6 (33) 
where the evolution matrix at NLO is given by 

U{n,nw) = U i0 \fi,fi W ) + ^ (jt> (0 WuO - U (0) (^fiw)j) (34) 

where is obtained by the 6x6 LO anomalous dimension matrix and J is obtained by 
the NLO anomalous dimension matrix. The explicit forms of these matrices can be found 
for example, in [42]. Since the O g contribution to B — > <pK s is order a s suppressed in 
the matrix element the Wilson coefficient C g (fx) should include, for consistency, only LO 
corrections: 

C g (fi) = U°(^fi w )C g (fi w ) (35) 

where U°(fj,, fj, w ) is obtained by the 8x8 anomalous dimension matrix of LO. 

The anomalous dimension matrix at NLO does depend on regularisation scheme. To 
avoid this problem, QCD corrections are carefully included in the literature [28]. As a 
result, the matrix element of B° — > <f)K° process is given by the effective Wilson coefficient, 

<^°|WSLiP°> (36) 

where 



G 



F , 



n eS= ~^ V tb v ts 



(37) 



.i=3 i=3 

The detailed expression of the effective Wilson coefficient can be found in [28]. The 
effective Wilson coefficient Cf (Cf ) includes all the QCD corrections mentioned above. 
We must emphases that these corrections also include the contribution from the chro- 
momagnetic type operator O g given in Eq. (32). Note that the LLO Wilson coefficient 

C< SM > itself is an order of magnitude larger than the others, however it enters as a QCD 
corrections so that a s /(47r) ~ 1/50 suppressed. As a result, the effect of O g in SM is less 

than 10% level in the each effective Wilson coefficients C^i^^- However, we will show 
that in supersymmetric theories, the Wilson coefficient for the operator O g {O g ) is very 



large and its influence to the effective Wilson coefficients Cf^g^^^ and Cf^^"^ are 
quite significant. 

Employing the naive factorisation approximation [43], where all the colour factor N 
is assumed to be 3, the amplitude can be expressed as: 



G 



The matrix element is given by: 



i=3 i=3 



{<t>K°\Oi\B Q ). (38) 



(<PK \O 3 \B ) = tx, (39) 

(<PK°\0 4 \B°) = ix, (40) 

{</>K°\0 5 \S°) = X, (41) 

(0K°|O 6 |5°) = l -X (42) 



with 



X = 2F^ K {ml)Um^p K • e,). (43) 

where Ff^ K (m^) is the B — K transition form factor and is the decay constant of 
meson. Note that the matrix elements for 0^ g ) and Oi( g ) are same for B — > (pK 
process. We use the following values for the parameters appearing in the above equation, 



m = 1.02 GeV, U = 0-233 GeV, (p K • e ) = (m 2 B - m\ + m 2 )] - m 2 ~ 13 

GeV, and F®^ K {mfy = 0.35 [44]. Finally, we discuss on the matrix elements of the 
chromomagnetic operator O g which is given by: 

imo 9 m = _^p> + 75 )^6^ (v^*) (44) 

where q^ is the momentum carried by the gluon in the penguin diagram. As discussed 
above, this contribution is already included in Cf£g, which in fact, is possible only when 
the matrix element of O g is written in terms of the matrix element of 03^6- It is achieved 
by using an assumption [28]: 

f = ^ (45) 

where (q 2 ) is an averaged value of q 2 . We treat (q 2 ) as an input parameter in a range of 
ml/A < (q 2 ) < ml/2. As we will see in the next section, our results are quite sensitive to 
the value of (q 2 ). 

4 Supersymmetric contributions to B — > (f)Ks decay 

As advocated above, the general amplitude A(<f>K) can be written as 

~A(<PK) = A SM ((f)K)+A°((f)K)+ ~A k± , (46) 
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Figure 1: The SM contribution (a) and the gluino-down squark contributions (b)-(f) to the 
B — > decay. 



where v4 SM , v4^, and A* refer to the SM, gluino, and chargino contributions, respectively. 
In our analysis, we consider only the gluino exchanges through AB = 1 penguin diagrams 

SUSY 

which give the dominant contribution to the amplitude A {<j>K). In Fig. 1, we exhibit 
the leading diagrams for B — > <pKs decay. At the first order in the mass insertion 
approximation, the gluino contributions to the Wilson coefficients G^ g at SUSY scale Ms 
are given by 



C 3 (M S ) 
C 4 (M S ) 
C 5 (M S ) 
C 6 (M S ) 
C 9 {M S ) 



(jFvtbVt s mq 

<^FVtbVts m a 



I 



-Bi(x) 



-Bo(x) 



1 



7 



I 



V2a 2 s 



r v v*™* K LLm 



-B x {x) + -B 2 (x) + -P x {x) + -P 2 (x 
^-B 1 (x) + -}-B 2 (x) - ^-Pi(x) - ]-P 2 {x 
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I 



(47) 



V2a 2 



2 ($LLm 



G F V tb V t *mq 
~G F V tb V t * s m 2 q 



-\ Bl {x) + 7 -B 2 {x) + \p 1 {x) A 2 {x) 

6 o o 2 



nig (1 



(6 d LL ) 23 ( -M 3 (x) + 3M 4 (x)) + (6 d LR h^ [-M,(.r) ■: 3 \U.V 



and the coefficients Cj j9 are obtained from Cj iS by exchanging L <-> R. The functions 
appear in these expressions can be found in Ref.[45] and x = rri~/m~. As in the case 
of the SM, the Wilson coefficients at low energy Cj i9 (/i), fi ~ 0(mb), are obtained from 
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Ci t9 (Ms) by using Eqs. (33) and (35). 

The absolute value of the mass insertions (Sjy B ) 23 , with A, B = (L,R) is constrained 
by the experimental results for the branching ratio of the B — > X s j decay [45,47]. These 
constraints are very weak on the LL and RR mass insertions and the only limits we have 
come from their definition, \(Sf LRR ) 23 \ < 1- The LR and RL mass insertions are more 
constrained and for instance with rrig ~ rriq ~ 500 GeV, we have \(Sf RRL ) 23 1 < 1.6 x 10~ 2 . 
Although the LR(RL) mass insertion is constrained so severely, as can be seen from the 
above expression of C g (Ms), it is enhanced by a large factor mg/rrib- We will show in 
the following that this enhancement plays an important role to reproduce the observed 
large deviation between sin 2/5 and S^k s - We should recall that in the supersymmetry 
analysis of the direct CP violation in the kaon system, the same kind of enhancement by 
a factor rrig/m s makes the LR and RL mass insertions natural candidates to explain the 
experimental results of e'/e [46]. 

As shown in Eq. (20), the deviation of S<pK s from sin 2(3 is governed by the size 
of |A SUSY |/]A SM |. Thus we start our analysis by discussing the gluino contribution to 
|A SUSY |/|A SM |. We choose the input parameters as 

m q = 500GeV, x = 1, q 2 = m 2 b /2, fj, = 2GeV (48) 

then, we obtained 

^SUSY 



0.13 (S d LL ) 23 + 55.4 (5 d LR ) 23 + 55.4 (5 RL ) 23 + 0.13 (5 d RR ) 23 . (49) 



The largest theoretical uncertainty comes from the choice of q 2 . We find that the smaller 
values of q 2 enhance the coefficients of each mass insertions and for the minimum value 
q 2 = ml/A gives 

4 SUSY 

-p^ ~ 0.23 (5 d LL ) 23 + 97.4 (5 d LR ) 23 + 97.4 (5 d RL ) 23 + 0.23 (5 d RR ) 23 . (50) 

Using the constraints for each mass insertions described above, we obtain the maximum 
contribution from the individual mass insertions by setting the remaining three mass 
insertions to be zero: 

I A SUSY I I a SUSY I 

H^<0.23, ^gpUl.56. (51) 

It is worth mentioning that (Sf R ) 23 and (S RL ) 23 contribute to S^x s with the same sign, 
unlike their contributions to e'/e. Therefore, in SUSY model with (Sf^^ ~ (S RL ) 23 , we 
will not have the usual problem of the severe cancellation between their contributions, but 
we will have a constrictive interference which enhances the SUSY contribution to S^Ks- 

Now let us investigate whether any one of the mass insertions can accommodate the 
observed large deviation between Sj/^k s and S^Kg- We start with LL (same for RR) 
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contribution. As can be seen from Eq. (20), a choice of the strong phase cos<5i 2 = ±1 
gives the largest deviation between S^Ks an d sin 2/5. Using the measured central value of 
sin 2(3, \(Sll)23\ — 1 an d cos#i 2 = ±1 in the CP violation master formula Eq. (18), we 
obtain S^k s m terms of 8iig(Sf jL / RR \)23 as: 

0.73 ± 0.45 sin(argfc (Rfl) ) 23 + 0.818) + 0.051 sin(2 argfe (RR) ) 23 + 0.818) 
^ 1.05 ± 0.45 cos(arg(^ L( ^ ) ) 23 ) 

(52) 

Then the minimum value of S^k s is obtained by sin(arg(^ L ( i?i? - ) ) 23 ) = =(=0.97 as 

S+Ks = 0.35. (53) 

We find that if the experimental value for S^Ks remains as small as the current central 
value, the SUSY models with LL or RR mass insertion can not provide an explanation 
for that. 

Next, we show that on the contrary, the (Sf jR / RL \)23 contribution can deviate S^Ks from 
sin 2/3 much more significantly. For \(S'[ R ^ RL ^)23\ = 0.01 and cos5i 2 = ±1, the mixing CP 
asymmetry is expressed as 

0.73 ± 1.95sin(argfe (RL) ) 23 + 0.818) + 0.95 sin(2 argfe (RL) ) 23 + 0.818) 
<pKs ~ 1.95 ± 1.95 cos(argfe (RL) ) 23 ) 

(54) 

and the minimum value is obtained with sin(arg(<5^ / j( i?L - ) ) 23 ) = (=f0.075, =f0.63) as 

S<f>K s — — 1- (55) 

In Fig. 2, we present plots for the phase of (Sll(rr))m an< ^ ^lr{rl))^ versus the mixing 
CP asymmetry S^,k s f° r cos5i 2 = I. We choose the three values of the magnitude of 
these mass insertions within the bounds from the experimental limits in particular, from 
B — > Xg'-f. Each plot shows a contribution from an individual mass insertion by setting 
the other three to be zero. As can be seen from these plots, the LR (same for RL) gives 
the largest contribution to S^Ks- In order to have a sizable effect from the LL or RR, 
the magnitude of (Sll(rr))23 nas ^° De °f or d er one and furthermore, the imaginary part 
needs to be as large as the real part. In any case, it is very difficult to give negative value 
of S^Ks from (5^ L ) 23 or (5^ i? ) 23 mass insertion. On the contrary, even if we reduce the 
magnitude of (5^)23 to the half of its maximum value, S<f,K s can still reach to a negative 
value. We also find that in the case of 1(5^/^)231 = 0.01, the minimum value of S<j,K s can 
be achieved without large imaginary part. Actually, we should notify that smaller values 
of the ratio x = m^/rn^ can lead to large LL and RR contributions to S^k s - However, the 
negative value of S^,k s can be only achieved by the light gluino mass around the current 
experimental limit and very heavy (order of TeV) squark mass. 

So far, we have only considered the cases where the strong phase 5u is given by 
cos5i 2 = ±1 so that the direct CP asymmetry C^k s was identically zero. However, since 
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ar g[(ii R(R I))23] 



Figure 2: The mixing CP asymmetry as function of arg[(<r^g)23] for three values of the | (^^4^)23 1 
where AB = LL(RR)(a), LR(RL)(b). The strong phase 5\2 is fixed at cos £12 = 1. 



the correlation between S^Ks and C^k s would be very useful to give some constrains on 
SUSY parameters, here we include arbitrary strong phases and demonstrate an example 
of SfyKs — C^k s correlation. By using the expanded formulae for S^Ks an d C^Ks in Eqs. 
(20) and (21), we find that for any value of #12, the plot of S$k s versus C^k s becomes an 
ellipse with its size proportional to smOsusV- 



(S^ Ks -sm2(3) 2 C, 



2 



cos 2 2(3 (2 sin 6 SUSY (2 sin 9 S usy 



= l / 56 ) 

/ I 4SUSY| \ 2 

\A SM \ 



In Fig. 3, we depict an example of the plot with |v4 SUSY |/|v4 SM | ~ 0.5 and Osusy — 
(7r/20, 7r/4, 7r/2, 37r/5, 47r/3). Since we used the full formulae in Eqs. (18) and (19) to 
create this figure, different Osusy does not give precisely rescaled ellipses. Nevertheless 
we can see the qualitative feature. The strong phase 5\2 = corresponds to the point at 
the right most tip of the ellipse. As £12 increases, it runs anti-clockwise and finishes a 
round when £ 12 = 2tt. Note that with |yl SUSY |/|y4 SM | ~ 0.5, a typical branching ratio is 
about 10 x 10~ 6 , which is within the range of experimental value. As experimental errors 
will be reduced, this kind of plot would provide very interesting constraints on SUSY 
parameters. 



5 CP asymmetry S^Ks m explicit SUSY models 

In this section, we study the CP asymmetry of the B — > <f)Ks in some specific SUSY 
models. We consider the minimal supersymmetric standard model (MSSM) (where mini- 
mal number of superfields is introduced and R parity is conserved) with the following soft 
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SUSY breaking terms 

V SB = + e ab (A^H b 2 qlu* Rj + A^H^d^ + A^H^e^. 

- BiiElE\ + H.c.) - -(m 3 gg + m 2 Ww a + m x BB) , (57) 

where i,j are family indices, a, b are SU(2) indices, e a b is the 2x2 fully antisymmetric 
tensor, with 6x2—1, and <j) a denotes all the scalar fields of the theory. We start with min- 
imal supergravity model and then we consider general SUSY models with non-universal 
soft breaking terms. We also discuss the impact of the type of Yukawa couplings on the 
prediction of the later model. 

5.1 minimal supergravity model 

In a minimal supergravity framework, the soft SUSY breaking parameters are universal 
at GUT scale, and we can write 

m 2 0a = m 2 , m t = m 1/2 , A% = A e^ A . (58) 

In this model, there are only two physical phases: (ft a = arg(^4*m!/ 2 ) and M = arg^m^)- 
In order to have EDM values below the experimental bounds, and without forcing the 
SUSY masses to be unnaturally heavy, the phases 4>a and M must be at most of order 
KT 1 and 1(T 2 , respectively [22]. 

It is clear that this class of models, where the SUSY phases are constrained to be very 
small and the Yukawa couplings are the main source of the flavour structure, can not 
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generate sizable contributions to the CP violating parameters. And it is indeed the case 
for S^Ks- We find that it is impossible to have a large deviation between S^Ks an d sin 2/5 
within the minimal supergravity framework. 

In fact, we find that even if we ignore the bounds from the EDM, and allow large 
values for SUSY phases, 4>a,h — 7r/2, still the SUSY contribution to S^, Ks is negligible. 
The suppression is mainly due to the universality assumption of the soft breaking terms. 
For instance, with ~ m ~ A ~ 200 GeV we find the following values of the relevant 
mass insertions: 

(5 d LL ) 2 3 ^ 0.009 + i 0.001, (59) 
(5 d RR ) 2 3 ~ -2.1 x 10~ 7 - i 2.5 x 10" 8 , (60) 
(5 d LR ) 2 3 ^ -2.5 x 10 5 — i 1.9 x 10^ 5 . (61) 

Clearly these values are much smaller than the corresponding values mentioned in the 
previous section and it gives negligible contributions to the CP asymmetry S^k s - Indeed, 
we find that the total S^k s m this example is given by S^k s = 0.729, which is about the 
value of sin 2/5. 



5.2 SUSY models with non— universal soft terms 

Now we consider SUSY models with non-universal soft terms. In particular, we focus 
on the models with non-universal A-terms in order to enhance the values of (Sf jR )23 and 
(S RL )23, which may give the dominant contributions to the CP asymmetry S,pK s - However, 
non-observation of EDMs leads to restrictive constraints on the non-degenerate A-terms 
and only certain patterns of flavour are allowed, such as the Yukawa and A-terms are 
Hermitian [23], or the A-terms are factorisable, i.e., (Y A )ij = AY or Y.A [46]. In the 
case of factorisation, the mass insertion (5 d R )n is suppressed by the ratio md/rriq. Here 
we will consider this scenario with the following trilinear structure: 



/ a 



A 



m 



a a 
b b b 

c c c 



(62) 



As pointed out in Ref. [24], in the case of non-universal soft breaking terms, the type of 
the Yukawa couplings (hierarchical or nearly democratic) plays an important role and has 
significant impact on the predictions of these models. If we consider the standard example 
of hierarchical quark Yukawa matrices, 



Y d 



— — diag (m u , m c , m t ) , 
v sin p 



v cos (5 



if. diag (m u , m c , m t ) .K + , 



(63) 
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where K is the CKM matrix, the relevant LR mass insertion is given by 

^ * ^ (* + y ^ *)- (64) 

where = Y^A d ^. It is clear that the dominant contribution to this mass insertion is 

given by the term ^22(^^)23-^33 which is still suppressed by the small entry of Y^. The 
non-universality of the squarks can enhance the LL and RR mass insertions, however this 
non-universality is severely constrained by the experimental measurements of AMk and 
Ek- Therefore, with the hierarchical Yukawa couplings we find that the typical values of 
the relevant mass insertions are at least two order of magnitude below the required values 
so that contributions to split the CP asymmetries Sj/^k s from S^Ks are again small. 

Now we consider the same SUSY model but with converting the above hierarchical 
Yukawa matrices to democratic ones, which can be obtained by a unitary transforma- 
tion. As emphasised in Ref . [24] that these new Yukawa textures (and their diagonalising 
matrices S^' d R ) have large mixing, which has important consequences in the SUSY re- 
sults. Thus, the element (YJ^) 23 has no suppression factor as before and the magnitude of 
(Sf R )23 can be of the desired order. As a numerical example, for m = = 200 GeV, 
(i.e., rriq ~ m g ~ 500 GeV) and assuming that \Aij\ G [m ,4mo] while the phases of the 
A-terms are chosen such that the bound of the EDMs are satisfied, we find that it is quite 
natural to obtain the following values of the mass insertion (6^)23- K^ljj) 23 1 — 0.005 and 
Arg[(^) 23 ] ~ 1.2 which leads to S 4>Ks ~ -0.2. 

6 Conclusions 

In this paper, we have studied the supersymmetric contributions to the CP asymmetry 
of B — > 4>Ks process. Using the mass insertion approximation method, we have derived 
model independent limits for the mixing CP asymmetry S^Ks- We found that the LR or 
RL mass insertion gives the largest contribution to S^Ks, while the LL or RR contribution 
is small. Thus, we conclude that if the deviation between S^k s an d Sj/^k s observed by the 
B-factory experiments (Belle and BaBar) remains as large as its present central value, 
the SUSY models with large (~ 10~ 3 ) LR(RL) mass insertions will be an interesting 
candidate to explain this phenomena. 

The Belle collaboration observed non-vanishing direct CP asymmetry C^k s which can 
be obtained only by simultaneous non-vanishing strong phase and SUSY CP violating 
phase. Thus, we studied the impact of the strong phase in our results for S^k s - We 
provided an example of a plot to show a correlation between S^Ks an d C$k s - O ur result 
will be useful to give constraints on some SUSY parameters as experimental errors will 
be reduced. 

We also applied our results to the minimal supergravity model and SUSY models with 
non-universal soft terms with two types of Yukawa couplings, hierarchal and nearly demo- 
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cratic Yukawa textures. We showed that only in SUSY models with large Yukawa mixing, 
SUSY contributions could be enhanced and reach the desired values to give significant 
contributions to the CP asymmetry S^Ks- This result motivates the interest in SUSY 
models with non-universal soft terms and also shed the light on the type of the Yukawa 
flavour structure. 

Note added in proof 

We have found an error in our program and the RR contribution is enhanced by one 
order of magnitude in this revised version. As a result, we agree with the results in Refs. 
[48,49]. 
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